Collision of two spin polarized fermionic clouds 
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We study the collision of two spin polarized Fermi clouds in a harmonic trap using a simulation 
of the Boltzmann equation. As observed in recent experiments we find three distinct regimes of 
behavior. For weak interactions the clouds pass through each other. If interactions are increased 
they approach each other exponentially and for strong interactions they bounce off each other several 
times. We thereby demonstrate that all these phenomena can be reproduced using a semiclassical 
collisional approach and that these changes in behavior are associated with an increasing collision 
rate. We then show that the oscillation of the clouds in the bounce regime is an example of an 
unusual case in quantum gases: a nonlinear coupling between collective modes, namely the spin 
dipole mode and the axial breathing mode which is enforced by collisions. We also determine the 
frequency of the bounce as a function of the final temperature of the equilibrated system. 

PACS numbers: 03.75.Ss, 03.75.Hh 



Spin transport has become in recent years a much stud- 
ied field in solid state systems, for example in mesoscopic 
phenomena, where the electronic spin degree of freedom 
is used to create new devices. Understanding the spin 
relaxation, diffusion and other transport properties is of 
fundamental importance in such fields. In atomic gases 
it has been studied mainly in spinor Bose gases [3, Q. 
There has been a renewed interest in spin transport in 
Fermigases which are a clearer parallel to electronic sys- 
tems |3H8[ . An important advantage of cold gases in such 
studies is the simplification due to the absence of relax- 
ation mechanisms for spin currents apart from direct col- 
lisions between atoms of different spin, unlike e.g. in a 
solid where collisions with the ionic lattice can be im- 
portant. In addition, the atomic interaction and initial 
temperature of the clouds are easily tuneable parameters. 
Finally, as we shall see, very large spin polarizations can 
be easily created leading to large spin currents. 

Here we study the collision of two clouds with oppo- 
site spin polarization following the recent experiments of 
A. Sommer et al. Q- We confine ourselves to the semi- 
classical regime, using a Boltzmann equation simulation. 
In contrast, a recent theoretical study has instead used a 
hydrodynamic approach based on a many-body equation 
of state [§] . 

One of the most striking experimental observations was 
the bouncing of the clouds off each other. We will here 
demonstrate that this phenomenon can be understood 
purely in terms of semiclassical collisions, without re- 
course to e.g. mean fields or other more complicated ef- 
fects. For instance, our approach predicts that all quan- 



tities considered here depend only on the square of the 
scattering length and not on its sign. Also, as we will 
show, the bouncing oscillations can be understood as an 
example of a nonlinear coupling between collective modes 
which, to our knowledge, has never before been studied 
in Fermi gases [Icj . 

We consider a system of two-component fermions la- 
belled by the spin index s = {y, ^} with equal mass m. 
We set h = ks = 1 throughout. Fermions of oppo- 
site spin can interact via s-wave collisions and the cross 
section is given by a = 4na 2 /(l + a 2 p;: cl /4), where a 
is the scattering length and p rc i the relative momen- 
tum of the two atoms. The fermions are confined in 
a cigar shaped harmonic trap with potential V(r) = 
|m(uja: 2 + ^yV 2 + ^z 2 ), where u> z < ui x = ui y . 

We assume that the system is in the normal phase 
and that the temperature is sufficiently high, so that the 
two spin distributions can be described semiclassically 
in terms of functions / s (r, p,t). The initial distribu- 
tions are created by sampling a Fermi-Dirac distribution 

f(t = 0) = ( e (p 2 /2m+V(r)- AI )/T lmt + norma li ze d t 

N/2 = N-f = N±, where Ti n ;t is the temperature and /i 
the chemical potential. The atom number N s for each 
species is held fixed during the simulation. Then each 
spin distribution is displaced in opposite directions along 
the z-axis by do/2. In the subsequent time evolution, 
the clouds begin to move towards the center under the 
harmonic trapping force resulting in a collision between 
them. After a sufficiently long time the center of mass 
energy will be transformed into the internal energy of the 
gas and a new equilibrium state will be reached, charac- 
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terized by a Fermi-Dirac distribution with temperature 
7fi na i (and a chemical potential /Ufinai)- Note that Tfi na i 
is a function of only the atom number, T; n ; t and c? an d 
can be calculated exactly from these values using energy 
conservation. 

Simulations were carried out for N = 10000, a range 
of initial cloud temperatures 0.2 < T^t/Tp < 10, inter- 
action strength 0.5 < |fcj?a| < 10, and initial distances 
between the centers of mass of the two clouds 0Aa z < 
do < I6cr z , where Tf = kp/2m = (SNujxUjyUJz) 1 ^ 3 and 
<j z = \J 'Tinit / 'muj2 . Our numerical setup is similar to 
the one described in ll|. Here we will only give a brief 



summary with a more detailed paper to follow 12j. 

The time-evolution of the distribution function 
/ s (r, p, t) is given by the Boltzmann equation, 

d t f s + (p/m) • V r / S - V r V ■ V p f s = -I[f 3 , fcl (1) 

where the left-hand side represents the propagation of the 
atoms in the potential and the right-hand side stands for 
the collision integral, 
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-f s m-h){i-h)] 



(2) 



The indices s and s label the two colliding atoms, the 
primed variables refer to quantities after the collision and 
O is the solid angle between the incoming and outgoing 
relative momenta. 

In order to simulate the Boltzmann equation we in- 
troduce a discrete time step. During each time step the 
atoms propagate following their classical trajectories. At 
the end of each time step collisions between the atoms 
are evaluated. The point-like atom picture is a discrete 
approximation of the continuous distribution function 
/ s (r,p,t). In order for this approximation to be accu- 
rate we represent each fermion by several test particles 
[TO . []]|]. The higher the ratio N/N of test particles to 
atoms, the more precisely the continuous distribution will 
be approximated. In this work we use N = 10N, which 
is sufficient for the range of parameters considered. 

The trajectories of the test particles are the same as 
the trajectories of the atoms and are given by the solu- 
tion to the classical harmonic oscillator equations. We 
evaluate collisions in the same way as described in [11|. 



First we test whether a pair of test particles fulfills the 
classical conditions for scattering. If this is the case we 
calculate the quantum mechanical scattering probability 
given by the Pauli term (1 — f' s )(l — fj) and accept or re- 
ject the collision accordingly. If a collision is accepted we 
calculate the new phase space coordinates of the particles 
respecting angular momentum conservation (l4| . 

We performed several tests of the simulation. We en- 
sured that the system thermalises to the correct equi- 
librium distribution (the Fermi-Dirac distribution in the 
presence of Pauli blocking, or the Maxwell-Boltzmann 



distribution if the quantum mechanical scattering prob- 
ability is set to one). We also checked the frequencies 
and the damping behavior of several collective modes, 
namely the dipole, the breathing and the quadrupole 
mode. A very important check is that the equilibrium 
collision rate measured during a simulation matches the 
theoretical prediction, as this will be crucial to obtain the 
correct damping behavior for the oscillation of two col- 
liding clouds. We performed tests both with and without 
Pauli blocking and obtained very good agreement with 
the values obtained analytically from the Fermi-Dirac 
distribution. 

The behavior of the clouds during the simulation can 
be studied by measuring the distance between their cen- 
ters of mass d(t) — (z^ — z{)(t). As in pj, we find 
that before they come to rest at thermal equilibrium, 
the motion of the clouds exhibits three typical behav- 
iors, see Fig. [T] and [l5|. Transmission: for sufficiently 
high temperatures and small interactions, the clouds os- 
cillate through each other (i.e. d(t) crosses zero at short 
times) with decreasing amplitude. Bounce: at low tem- 
peratures and strong interactions the clouds bounce off 
each other several times (in each bounce the motion of 
the center of mass of each cloud is reversed at short times 
and without d(t) crossing zero) before a longer period of 
slow approach. Intermediate: between the transmission 
and bounce regimes there is a range of temperatures and 
interactions where the slow approach behavior is visible 
from the start and neither bounces nor transmissions are 
observed. 

The dependence of the behavior on temperature and 
interactions is related to the variation in the collision rate 
7 in the overlap region between the two clouds. As the 
collision rate decreases the system behavior changes from 
the bounce regime to intermediate, and finally to the 
transmission regime. From the top panel of Fig. Q] we see 
that in the bounce regime the oscillations in the collision 
rate integrated over a volume in the overlap region follow 
closely the oscillations of d(t), whereas no such variation 
is apparent in the transmission regime. In addition, we 
can compare the collision rate with the typical timescale 
for macroscopic motion (uj^ 1 ). We see that the gas is 
strongly hydrodynamic in the bounce regime (uj z /j *C 
1) and becomes collisionless in the transmission regime 

(w,/7~l). 

Figure [2] shows the transitions between these regimes 
in the (T/Tf) final i I^Ffllfinai plane for a fixed aspect ratio 
and different values of do in units l z = 1/ \/muj z . Here 
the Fermi temperature is defined as Tp — k F /2m — 
(37r 2 no) 2 / 3 /2m, where no is the total atomic density 
no = n^o + n^o in the trap center. Note that as the den- 
sity is a function of the temperature, Tp and kp change 
during the simulation and hence kpa varies during the 
evolution. The quantities given here are equilibrium val- 
ues for t — > oo. 

The transition between transmission and intermediate 



FIG. 1. Bottom panel: the normalized dipole mode d(t)/do (solid lines) and breathing mode b(t)/boo (dashed lines) for the 
three different behaviors: transmission (left), intermediate (middle) and bounce (right). Top panel: the corresponding collision 
rate per particle y/uiz measured in the region with \z\ < a z /2 around the trap center, where a z = \jT\^itjmLo%. 
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FIG. 2. (Color online) The transition between bounce and 
intermediate regimes (filled symbols, lines to guide the eye) 
and between intermediate and transmission regimes (empty 
symbols). Red circles correspond to do = 43. 1Z Z , blue trian- 
gles to do = 64.6Z Z and green squares to do = 129. 3l z , where 
l z — 1/ yfmuj z . It is clearly visible that the intermediate- 
transmission transition is independent of do. The dashed line 
corresponds to constant relaxation time l/rdi p = 1.83oj z . 



regimes is defined as d(t) reaching but not crossing zero 
at short times. It is clearly visible from Fig. [5] that this 
transition is independent of do and hence can be under- 
stood entirely from the final equilibrium properties of the 
system. More precisely, it can be understood as a con- 
sequence of the change in the relaxation time Tdi P of the 
spin dipole mode of the do — system, which is closely re- 
lated to the collision rate per atom. For sufficiently high 
temperatures T > Tp the relaxation time can be calcu- 
lated from the Maxwell-Boltzmann distribution [IE| and 
equals 
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(3) 



where f(y) = 1 — y + y 2 e y T(0,y). For sufficiently low 
collision rate (w z Tdi P ~> 1), the gas can be said to be 
collisionless and therefore the clouds undergo indepen- 
dent oscillations without interacting strongly with each 
other. For l/raip > 2w 2 the dipole mode is overdamped 
IE El- To compare with the simulations, we calculate 
the value of Tdi p for the various points lying on the curve 
separating the transmission and intermediate regimes of 
Fig- HI using Eq. (|3]). We find that they lie on the curve 
of constant l/r^p = 1.83(l)u; z . 

The transition between the intermediate and the 
bounce regime is defined to occur when the first bounce 
ceases to reverse the motion of the clouds, or in other 
words when d(t) ceases to have a minimum and becomes 
a monotonically decreasing function of t. This transition 
depends on do- In the bounce regime we typically see 
an initial strong collision followed by oscillations of d(t) 
which eventually die out as d(t) —> 0. This oscillatory 
behavior continues into the intermediate regime. The 
bottom panel of Fig. Q] shows plots in the three regimes 
of the normalized amplitudes of the dipole mode d(t)/do 
and the breathing mode &(£)/6oo where b(t) = (zl+zf)(t). 
In the bounce regime, the frequency of d(t) is identical to 
the frequency of b(t), see Fig. Q] (right), and suggests the 
existence of a nonlinear coupling between the two modes. 
However, in the transmission regime, see Fig. [T] (left), the 
oscillation frequency of d(t) becomes closer to that of the 
dipole mode of the non-interacting gas w z , and so the 
spin dipole mode decouples from the breathing mode. 

For a more quantitative analysis we fit to the function 



d(t) = Be- 1 ' ( 1 + De- t/E sin (wt 



(4) 



The first term is related to the spin-drag coefficient mea- 
sured in 0]. It dominates the overdamped behavior of 
d(t) at long times with a characteristic timescale C which 
we will analyze elsewhere 12| . The second term describes 
the coupling between the spin dipole and the breathing 
mode. 
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FIG. 3. (Color online) The frequency uj/ui z of the dipole mode 
d(t) (red circles) and the breathing mode b(t) (blue triangles) 
versus the final temperature for \k,Fa\ = 1. All data were ob- 
tained for equal initial temperature Ti n it = 0.4T> by varying 
do. The solid line is the prediction from [181 ]. 



The oscillations of the axial size of the cloud are fitted 
using the dependence 



b(t) /b QO = l + D'e' t/E ' sin (u't + </>'). 



(5) 



From the fit to our simulations, we observe that D 1 ~ 
2D, E' ~ E, uj' ~ u and </>' ~ </>. We interpret this 
result by describing the density profile of the gas by the 
ansatz n s (z,t) — an s (az ± /3,0), where a — a(t) and 
j3 — /3(t) represent the breathing and spin dipole modes 
respectively. Using this expression, we can calculate b(t) 
and d{t) as functions of a and (3 and in the weakly non- 
linear regime, we obtain d(t) ~ 2/3/a and b(t) ~ 1/a 2 . 
Assuming that j3(t) is overdamped and a(t) is a damped 
oscillator we then obtain the time dependence in Eqs. Q 
and (|5|) and the corresponding relations between the fit 
coefficients. 

We are also able to study the temperature dependence 
of the frequency of the two modes, see Fig. [3] We ob- 
tain the frequency by fitting to the functions (|4]) for the 
dipole mode and ([5]) for the breathing mode. Since the 
frequency shows a weak time dependence we ignore early 
times by imposing a cut-off d c — 2 • 0.7a z on the ampli- 
tude and fitting only times t > t c for which d(t) < d c . 
Since the function d(t) is not monotonic the correspond- 
ing time cut-off t c is not continuous with changing do, 
which leads to a small systematic error. We estimate 
this error from the cut-off dependence of w. The statisti- 
cal error of the fit is several orders of magnitude smaller 
and hence negligible. 

At low temperature, the common frequency of the spin 
dipole and breathing modes is close to the hydrodynamic 
prediction y/\2j5uj z « 1.55w z [3, Hit ], as observed ex- 
perimentally [7j. As we increase temperature, the fre- 
quency approaches 2uj z , the non- interacting value. At 



higher temperatures the spin dipole mode frequency be- 
comes ill-defined due to large damping. If we continue to 
increase the temperature, the damping becomes progres- 
sively smaller and the dipole mode frequency approaches 
lu z , the value for the non- interacting gas. We also com- 
pare our results with an earlier prediction for the fre- 
quency of the breathing mode |18[ . The two dependen- 
cies are close to each other although some discrepancy 
remains. We attribute it to the fact that to estimate the 
collisional rate in the cloud Ref. [HI neglects the Pauli 
principle and assumes a Gaussian phase-space density. 

In conclusion, we have studied the collision of two spin 
polarized fermionic clouds using a Boltzmann equation 
simulation. We found various regimes of behavior, char- 
acterized the transitions between them as a function of 
interaction strength and temperature, and related them 
to the collision rate in the overlap region between the 
clouds. In particular we showed that the bounces can 
be explained purely as a semiclassical collisional phe- 
nomenon, without the need for more complex many-body 
effects. We also demonstrated that the bounces are a rare 
example of nonlinear mode coupling, in which the spin 
dipole and breathing modes interact. In future we aim 
to extend this study to other closely related problems 
such as the collision of clouds with unequal populations 
[20j j . between clouds of atoms with unequal masses plj ]. 
and to the Fermi liquid regime at T <C Tp, using the 
Landau-Boltzmann equation. 
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